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Abstract 

Using first-principles calculations we show that the band gap of bilayer sheets of semiconducting 
transition metal dichalcogenides (TMDs) can be reduced smoothly by applying vertical compressive 
pressure. These materials undergo a universal reversible semiconductor to metal (S-M) transition at 
"^ , a critical pressure. S-M transition is attributed to lifting the degeneracy of the bands at fermi level 

caused by inter-layer interactions via charge transfer from metal to chalcogens. The S-M transition 
can be reproduced even after incorporating the band gap corrections using hybrid functionals and 
GW method. The ability to tune the band gap of TMDs in a controlled fashion over a wide range 
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I. INTRODUCTION 

Despite being the most promising two-dimensional (2D) material gapless graphene has 
limitations for its applications in nanoelectronics and nanophotonics. This led to the finding 
of other 2D-materials with finite band gap such as BN, transition metal dichalcogenides 
(TMDs) and transition metal oxides (MO2). BN sheet is an insulator and modifying their 
band gap for optical and electronic applications is still a challenge. Depending on the 
combination of metal and chalcogens, TMDs offer a wide range of 2D materials: metalsi*^, 
superconductors^!^, charge density wave systems^i^, Mott insulators^ and semiconductors^*^. 
Semiconducting two-dimensional TMDs include M0S2, MoSe2, MoTe2, WS2 and WSe2 and 
have emerged as promising materials for a range of applications. 

Few layers to monolayer of M0S2 and other TMDs have been successfully synthesized^"— 
and their optical absorption and photoconductivity have been studied^i^S. A single layer 
of TMDs having stoichiometry of MX2, consists of a hexagonally arranged transition metal 
(M = Ti, Nb, Ta, Mo, and W) sandwiched between two layers of chalcogen atoms (X = S, 
Se, and Te). Within a layer the metals and chalcogens form strong ionic-covalent bonds, 
whereas in bulk TMDs, the layers are bonded by weaker van der Waals (vdW) interaction. 
The bulk TMDs are indirect band gap semiconductors having band gaps in the range of 1.0- 
1.35 eV^. With the reduction in number of layers the band gap increases^ and transforms 
into a direct gap for a single layer TMD^"— . In order to use these materials as building 
blocks in nanoelectronics, their electronic properties need to be modified. This has been 
achieved by doping^i^ and intercalation^^"^. The TMD based field effect transistors with 
high room-temperature current on/off ratios^- and higher on- current density^^^^ as well as 
integrated circuits— have been successfully fabricated. 

Tuning band gap of 2D materials — for their potential application in electromechanical 
devices, tunable photodetectors and lasers — has been a challenge in band gap engineering. 
Applied strain or electric field offers a novel way of modifying the band gaps over a wide 
range. Even for graphene, it has been shown using density functional theory calculations 
that uniaxial strain on monolayers^ and applied vertical electric field to bilayers^"— open a 
small band gap. Similar, theoretical studies for TMDs show the semiconductor to metal (S- 
M) transition for bilayers subjected to vertical electric field (0.2-0.3 V/A)-^ and for mono- 
and bilayers under biaxial straiii^^. These techniques are promising but suffer from the 
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FIG. 1. (Color online) Structure of AB-stacked MX2 bilayer with M and X atoms are shown by 
purple and yellow spheres, respectively. Side view of the bilayers with AA and AB stacking. The 
blue dotted lines show the unit cells. 

practical applicability. For example, the electric field required to achieve S-M transition 
is too large, and a reversible way of applying a biaxial strain in 2D materials has yet to 
be demonstrated in laboratory. On the other hand it has been experimentally shown that 
the band structure of bulk TMDs can be modified by application of compressive strairt^. 
Here, we investigate the effect of normal compressive strain (NCS) on electronic properties 
of semiconducting bilayer TMDs: M0S2, MoSe2, MoTe2, WS2 and WSe2. The band gaps 
of these materials decrease gradually with the applied NCS. A reversible semiconductor to 
metal transition was observed after a threshold pressure Pth was achieved. The Pth depends 
upon the material as well as the stacking pattern of the two layers. The hybrid functional 
and GW methods were used to correct the PBE gap. The S-M transition was found to be 
independent of the methods. This offers a wide range ( 1.9 - 0.0 eV) of reversible band gap 
tuning, which can be utilized for various applications. 



II. COMPUTATIONAL DETAILS 



The calculations were performed using ab-initio density functional theory (DFT) in con- 
junction with all-electron projector augmented wave potentials'"^ and the Perdew-Burke- 
Ernzerhof' generalized gradient approximation (GGA) to the electronic exchange and corre- 



lation, as implemented in the Vienna Ab initio Simulation Package (VASP)^. We optimize 
the structure of bilayer TMDs (MX2 with M = Mo, W and X = S, Se, Te), using the unit 
cells as shown in Fig. [TJ A well converged Monkhorst-Pack k-point set (15 x 15 x 1) was used 
for the calculation and conjugate gradient scheme was employed to optimize the geometries 
until the forces on every atom were < 0.005 eV/A. Sufficient vacuum was used along z di- 
rection i.e. perpendicular to the 2D sheet, to avoid spurious interaction among the periodic 
images. The lattice parameter and interlayer distance of the optimized structures are listed 
in Table |T] and are in good agreement with the previously reported PBE values^. 

Structural Parameters (A) Band Gap (eV) 

Material Functional a c Bilayer separation Bulk Monolayer Bilayer 

AB AA AB AA 

PBE 3.18 12.42" 4.14 4.54 0.94 1.68 1.59 1.63 

M0S2 PBE + vdW 3.21 12.42 3.11 3.72 0.85 - 1.25 1.47 

Exp. S.ie^'O 12.29^ - - 1.23^ 1.9^ 1.6^ - 

PBE 3.32 13.06** 3.88 4.36 0.87 1.43 1.42 1.41 

MoSe2 PBE + vdW 3.33 13.06 3.19 3.78 0.87 - 1.18 1.41 

Exp. 3.30''0 12.94^ - - 1.09^ 

PBE 3.56 13.95** 3.95 4.97 0.67 1.06 1.03 1.05 

MoTe2 PBE + vdW 3.54 13.95 3.37 4.07 0.68 - 0.88 1.02 

Exp. 3.52^ 13.97-^ - - 1.00^° 

PBE 3.18 12.99** 4.31 4.56 1.25 1.81 1.77 1.79 

WS2 PBE + vdW 3.19 12.99 3.39 3.95 1.24 - 1.57 1.70 



Exp. 3.15^^2 12.3242 . . 1.358 



PBE 3.32 13.38** 4.11 4.43 1.09 1.53 1.51 1.52 

WSe2 PBE + vdW 3.34 13.38 3.35 3.98 1.08 - 1.43 1.51 



Exp. 3.28^2 12.9642 _ _ 1 208 



**Same values as vdW 

TABLE I. Structural parameters and band gaps of bulk, monolayer and bilayer of TMDs calculated 
using PBE, PBE+vdW. For comparison purpose the corresponding experimental values are also 
included. 



The week van der Waals interaction between the layers has an effect in determining 
the interlayer distance for the bilayers as well as for the bulk MX2. The van der Waals 
interaction originates from dynamical correlations between fluctuating charge distributions 
and cannot be described by the PBE functional. Consequently, the relaxed bilayer distance 
obtained by PBE is off by approximately 1 A compared to the bulk interlayer distance. 
We incorporate the van der Waals interactions through a pair-wise force field following 
Grimme's DFT-D2 method^. The parameters used in DFT-D2 method are thoroughly 
optimized for several of the DFT functionals, including PBE. We performed a comprehensive 
test on the reliability of the empirical parameter by calculating the bulk phase of TMDs. 
The lattice parameters obtained using this approach are in very good agreement with the 
experimental values^^^, Table [H The calculated DFT-D2 interlayer distances in bilayer 
TMDs have decreased significantly compared with PBE, Table [B 

III. RESULTS AND DISCUSSION 

The two layers in a bilayer TMDs, can be arranged in AA and AB stacking. In the AA 
(AB) stacking the Mo atoms in one layer is on top of the Mo (S) atoms in another layer. 
Fig. [H AB-stacking is preferred in both bulk and bilayer TMDs. In the case of a bilayer 
TMD the energy differences between two types of stackings are very small. Therefore, in the 
current work, the effect of applied normal compressive strain on the band structure of both 
AA- and AB- stacked TMDs, was studied. The strain was calculated as e = {dQ—d)/dQ, where 
do and d are equilibrium and instantaneous interlayer distance, respectively. A constrained 
relaxation scheme, where the positions of metal atoms remained unchanged, was used for 
the optimization of the structures. This prevented the relaxation of instantaneous structures 
to the original positions. All the reported results here onwards include the vdW interactions 
via DFT-D2 method, unless mentioned otherwise. 

The overall behavior of the band structure under the applied pressure remains similar for 
the both type of stackings. Therefore, here we will present mainly the results of AB-stacked 
TMDs. The PBE+vdW band structure calculations show that the unstrained bilayers of 
all these AB stacked TMDs are indirect band gap semiconductors. In comparison to mono- 
layer, the band gap of bilayer decreases slightly and each band becomes doubly degenerated 
as shown in Fig. [2] for M0S2. This clearly demonstrates the absence of chemical interac- 
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FIG. 2. (Color online) Band structure, total DOS, PDOS of Mo and S for M0S2 bilayer at (a) 0, 
(b) 0.18 and (c) 0.34 strain. 



tion between the layers. As the interlayer separation decreases, the layers start to interact 
chemically, which leads to lifting of double degeneracy of the bands. Fig. |2l This splitting 
increases with the increasing NCS. The valence band maxima (VBM) moves away/towards 
the fermi level at K-/ F-point. The S-M transition occurs when VBM crosses the fermi level 



at the r point after a critical applied NCS. Likewise, the conduction band minima (CBM) 
also moves towards the fermi level with increasing NCS. The band structures of the other 
MX2 materials undergo similar changes (lifting of double degeneracy, and S-M transition) 
under applied NCS. The interlayer distance at which the S-M transition takes place is listed 
in Table Hi] for the bilayers. 

In order to determine the constituents of the electronic bands, the total and projected 
density of states (PDOS) were calculated for the bilayers under various NCS and shown for 
M0S2 in Fig. m The total DOS shows a gap near the fermi level, which reduces due to the 
shift of the CBM and the VBM towards the fermi level with increase in strain. For the 
unstrained M0S2, the CBM and the VBM are constituted by Mo-d and S-p orbitals, which 
can be clearly seen from PDOS. With the increasing NCS, the contributions to CBM from 
Mo-d and S-p orbitals decrease, while from S-d increase. The change in the band gap as a 
function of applied NCS follows similar pattern for all the materials studied here, (Fig. [3]^a)). 
This change is reversible, i.e in the absence of the applied pressure the structure relaxes back 
to the original structure with the complete recovery of band gap, which is very important 
for sensor applications. Similar changes in band gap^ and resistances^ were reported for 
bulk M0S2 under applied pressure, however, unlike bilayers no S-M transition was observed 
even at very high pressure (40 GPa). 

In order to access the feasibility of S-M transition in experiments, the applied pressure 
(P) was calculated from the energy cost per unit area in reducing the interlayer distance by 

Material Bilayer separation at transition (A) Transition Pressure (Pth) (GPa) 
AB AA AB AA 



M0S2 


2.14 


2.60 


8.52, 


14.47 (HSE) 


8.37 


MoSe2 


2.23 


2.64 




8.37 


9.54 


MoTe2 


2.69 


2.93 




5.10 


8.71 


WS2 


1.80 


2.47 




16.28 


10.54 


WSe2 


2.21 


2.79 




15.83 


12.28 



TABLE II. Interlayer distance and pressure required for semiconductor to metal transition, Pth of 
TMDs for both AA and AB stacking calculated using vdW. 
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FIG. 3. (Color online) (a) Change in band gap with applied strain for MX2 bilayers for both AA 
and AB stacking using PBE+vdW method. The inset show the same using PBE+vdW, HSE and 
GW method for AB stacked M0S2 bilayer. (b) Average Bader charge of M and X for all the bilayers 
as a function of strain. The left and right panels corresponds to M = Mo and W, respectively. 
Same color and symbols are used for the respective bilayers as defined in (a). Isosur faces of charge 
accumulation (green) and depletion (red) of M0S2 bilayer under strain of (c) 0.18, (d) 0.24 and (e) 
0.31. 

Ad = {do — d) as per the following equation, 

(do -d)A ^ ^ 

where A is the area of the unit cell, E and Eq are the energies and d and do are the 
interlayer distances of the strained and unstrained bilayer. The pressure required for S-M 
transition, Pj/j for each material is listed in Table [TTl The calculated pressure range is easily 
achievable experimentally, hence makes tuning of band gap via NCS very attractive for 
various applications. Percentage reduction in band gap with increase in pressure is plotted 
in Fig. lU The plot is linear for most of the bilayers except for AB stacked WX2. For 
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MoTe2, the change in band gap is very less initially but increases drastically after a pressure 
of 1.3 GPa. The slope and hence the response of change in band gap to applied pressure 
for MoTe2 (AB) is the largest among all the TMDs. Surprisingly, S-M transition occurs at 
relatively lower pressure for the AA stacking than the AB except for MoTe2. For a given 
metal atom, the S-M transition pressure decreases as we go down the column of X atom 
in the periodic table. This is caused by increased delocalization of atomic orbitals, which 
leads to reduced interaction between M and X atoms resulting in S-M transition at lower 
pressure. Such behavior is consistent with the trend of band gap, which also decreases from 
StoTe, Table H 

Similar evidence also comes from the Bader charge analysis^"—. The average charge on 
M and X as a function of applied strain is shown in Fig. |3]^b). The charge transfer from 
M to X increases as a function of NCS, which is more prominent for WX2. Furthermore, 
for a given M atom, the value of average charge on X atom decreases as we go down the 
periodic table i.e., from S to Te, replicating closely the trend of electron affinity of X atoms. 
This, essentially, results into weaker interaction between M-X atom, which leads to the S-M 
transitions pressure to be the lowest for the Te and the highest for the S. Furthermore, in 
comparison to W, Mo has lower ionization potential, which measures the ability to donate the 
charge. Mo can donate charges relatively, easily and thereby facilitate the X-X interaction, 
which results in S-M transition in M0X2 at lower pressure compared with the WX2. 

These conclusions were further supported by redistribution of charges, which was calcu- 
lated by taking the difference of the total charge of bilayer and two isolated layers for M0S2, 
as shown in Fig. [3] c-e. The maximum redistribution of the charge occurs around the inner 
S atoms and the adjacent Mo-S bonds. As expected, for low strains (below e = 0.21) no 
charge rearrangement was observed. At the non zero strains, the charge started accumulat- 
ing on the inner S atoms as well as on adjacent Mo-S bonds. Fig. [3] c-e. However, the charge 
depletion was observed predominantly from the inner S atoms. With the increasing NCS, 
the amount of charge redistribution also increases, indicating the enhanced S-S interaction, 
which essentially causes the S-M transition. Due to smaller S-S distances the interaction is 
better for the AA-stacking compared with the AB. This is why the Pth for AA-stacking is 
lowered compared with the AB-stacking. 

In order to gain further insight, the band decomposed electron densities of the VBM 
and CBM, at different strains are plotted and shown in Fig. [51 Consistent with the above 
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FIG. 4. (Color online) Percentage change in band gap with applied pressure for MX2 bilayers 
obtained using PBE-vdW. Inset: The same plot using PBE-vdW and HSE method for AB stacked 
M0S2 bilayer. 

analysis, the VBM originates from the interaction of inner S atoms and adjacent Mo-S 
bonds, Fig. [5] a. With the increasing interaction caused by apphed NCS, the contribution 
from inner S atoms further increases. The VBM is extended and hence the effective mass 
does not change much even after the apphcation of NCS, implying maintenance of quality 
of conductivity. Furthermore the lowering of the symmetry of VBM is also evident. The 
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CBM completely originates from the Mo-S bonds, Fig. [5]b. Once again the nature is more 
extended and hence will keep the conductivities intact. The modification of gap without 
significant alteration of the dispersion of gap, is very importance for electronic applications. 

Due to the presence of artificial self-interaction--^ and the absence of the derivative dis- 
continuity in the exchange-correlation potential, DFT in the local-density (LDA) and gen- 
eralized gradient approximation (GGA), suffers from the underestimation of the band gap. 
The calculated PBE/GGA band gap of bulk, mono- and bi- layer M0S2 is underestimated 
by an amount of 0.29 eV(23.7%), 0.22 eV(11.5%) and 0.01 eV(0.5%) than their experimental 
values, respectively (Table [I]). The excellent agreement for the bi-layer is misleading as it 
corresponds to a wrong geometry (extremely large c value) obtained by PBE in the absence 
of vdW interactions. Band gap of a bilayer with the geometry optimized by PBE+vdW is 
underestimated by 0.33 eV(20.6%) with respect to the experimental value. For the com- 
parison, band gap using LDA functional was also calculated for M0S2. The LDA predicts 
correct band gap for monolayer but gives values 0.47 eV(38.4%) and 0.66 eV(4Ll%) lower 
than the experimental ones for the bulk and bilayer M0S2 respectively. 

In order to correct the GGA/LDA band-gaps, we used hybrid Heyd-Scuseria-Ernzerhof 
(HSE)^i^ functional. In the HSE approach, the exchange potential is separated into a long- 
range and a short- range part. The 1/4 of the PBE exchange is replaced by the Hartree-Fock 
(HF) exact exchange and the full PBE correlation energy is added. HSE is shown to correct 
the GGA band gaps^"^, significantly by partially correcting the self interaction. However, 
the results are system dependent. The calculated band gaps for M0S2 bilayer using HSE is 
1.64 eV which is in very good agreement with the experiment (1.60 eV). 

We also corrected the PBE band gaps by using DFT with many-body perturbation theory 
in the GW approximation^^. A partially self consistent GWO method was used in which the 
G was iterated but the W was kept fixed to the initial DFT (PBE) WO. A default cutoff 
(280.0 eV) for the wave function was used for the GW calculation. Convergence study was 
performed for number of bands, k-points and frequency grid points to achieve a convergence 
within 10-20 meV for the band gaps. A F-centered Monkhorst-Pack k-point grid of 18 x 18 x 1 
mesh was used for calculating GW band structure. Like HSE, the GWO also improves PBE 
band gap and overestimates it by 0.18 eV(11.3%) (for bilayer) compared with experimental 
value. The GWO band gap is slightly larger than the HSE. There are not many experimental 
results for the band gaps of the 2D-materials. Furthermore, the quality of as grown 2D- 
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sheets can also play an important role in band-gap determination. In order to get a more 
reliable comparison more experiments are required. The band gap calculated using HSE 
functional is the closest to the experimental value among all the methods used in this work 
for bilayer M0S2. Both HSE and GW calculations are computationally very expensive and 
hence we applied these methods only to M0S2. 

We checked the robustness of the S-M transition by calculating the band-structure of 
M0S2 as a function of normal compressive strain using HSE functional and GW method. 
The plot of HSE and GW band gap as a function of applied strain is compared with the 
one for PBE+vdW in Fig. |3] a inset. A comparison of percentage reduction in band gap 
with increase in pressure for different materials and methods is shown in Fig. |H Although, 
the calculated HSE and GW band gap of M0S2 was slightly larger than the PBE, but 
the nature of the band structure as well as the S-M transition remained unchanged. GW 
method does not give any force or energy of the systems and therefore pressure could not 
be calculated and compared with other methods for it. Since the band gap obtained from 
HSE and GW calculation was slightly more than the ones with the GGA, an increase in 
transition interlayer distance (as well as Pth for HSE) was observed for the S-M transition. 
The excellent agreement between the HSE and experimental band-gap, suggests that the 
transition temperature would be closer the HSE values. However, more experiments are 
required to claim a good numerical accuracy. 

IV. CONCLUSION 

In conclusion, we report reversible band gap engineering of semiconducting bilayer of 
transition metal dichalcogenides by applied normal compressive strain. The band gap can 
be tuned in a large energy range without modifying the conductivity significantly. Fur- 
thermore, a universal reversible semiconductor to metal transition was observed for all the 
semiconducting TMDs. The reduction of gap as a function of applied pressure is caused 
by interlayer interaction, which eventually lifts the degeneracies of bands and moves them 
closer to the fermi level. For a given M in MX2, the threshold pressure Pth, to achieve 
S-M transition decreases as X changes from S to Te. Furthermore, M0X2 has lower Pth, 
compared with the WX2. The PBE band gap of M0S2 was corrected using hybrid HSE 
functional and GW method. The S-M transition can be reproduced even after applying the 
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FIG. 5. (Color online) Isosurfaces of band decomposed charge density of (a) valence band maxima 
and (b) conduction band minima of AB staked M0S2 bilayer at applied strain of 0, 0.18 and 0.31, 
respectively. 

band gap corrections, however, the Pth changes to a larger value. The tantalizing possibility 
of reversible tuning of more than 1.9 eV of energy gap by applying compressive strain as 
shown in the present work, would make the usage of TMDs in wide range of applications 
spanning from sensors to electronics. 



* labhishek@mrc. iisc.ernet.ini 

1 Y. Ding, Y. Wang, J. Ni, L. Shi, S. Shi, and W. Tang, [Physica B, 406, 2254 (20lT) 



13 



2 A. Ayari, E. Cobas, O. Ogundadegbe, and M. S. Fuhrer, |J. Appl. Phys., 101, 014507 (2007)| . 
^ A. M. Gabovich, A. I. Voitenko, J. F. Annett, and M. Ausloos, 



Supercond. Sci. Tech., 14, Rl (2001) 



E. Morosan, H. W. Zandbergen, B. S. Dennis, J. W. G. Bos, Y. Onose, T. Klimczuk, A. P. 



Ramirez, N. P. Ong, and R. J. Cava, Nat. Phys., 2, 5 44 (2006), 
^ K. Rossnagel, O. Seifarth, L. Kipp, M. Skibowski, D. VoB, P. Kriiger, A. Mazur, and J. Poll- 



mann, 



Phys. Rev. B, 64, 235119 (2001) 



6 W. Z. Hu, G. Li, J. Yan, H. H. Wen, G. Wu, X. H. Chen, and N. L. Wang, |Phys. Rev. BTTG; 
045103 (20071 . 



■'' A. F. Kusmartseva, A. Akrap, H. Berger, L. Forro, and E. Tutis,' Nat. Mater., 7, 960 (2008) 
8 K. K. Kam and B. A. Parkinson, [J. Phys. Chem., 86, 463 (1982) . 



9 



S. W. Han, H. Kwon, S. K. Kim, S. Ryu, W. S. Yun, D. H. Kim, J. H. Hwang, J.-S. Kang, 



J. Baik, H. J. Shin, and S. C. Hong, phys. Rev. B, 84, 045409 (2011) 
10 R. F. Frindt and A. D. YofFe, |Proc. R. Soc. Lon. Ser.-A, 273^ 69 (1963)| 



11 



K. S. Novoselov, D. Jiang, F. Schedin, T. J. Booth, V. V. Khotkevich, S. V. Morozov, and 



A. K. Geim, |Proc. Natl. Acad. Sci. USA, 102, 10451 (2005) . 

12 R. F. Frindt, |J. Appl. Phys., 37, 1928 (1966) . 

13 P. Joensen, R. F. Frindt, and S. R. Morrison, pVIater. Res. Bull., 21, 457 (1986) . 

1^ A. Schumacher, L. Scandella, N. Kruse, and R. Prins, purf. Sci., 289, L595 (1993) . 

1^ J. N. Coleman, M. Lotya, A. ONeill, S. D. Bergin, P. J. King, U. Khan, K. Young, A. Gaucher, 
S. De, R. J. Smith, I. V. Shvets, S. K. Arora, G. Stanton, H.-Y. Kim, K. Lee, G. T. Kim, G. S. 
Duesberg, T. Hallam, J. J. Boland, J. J. Wang, J. F. Donegan, J. C. Grunlan, G. Moriarty, 
A. Shmeliov, R. J. NichoUs, J. M. Perkins, E. M. Grieveson, K. Theuwissen, D. W. McComb, 



P. D. Nellist, and V. Nicolosi, [Science, 331, 568 (2011) . 
16 B. K. Miremadi and S. R. Morrison, |J. Appl. Phys., 63, 4970 (1988) . 
i'^ D. Yang and R. Frindt, |J. Phys. Chem. Solids, 57, 1113 (1996) . 
1^ H. S. S. Ramakrishna Matte, A. Gomathi, A. K. Manna, D. J. Late, R. Datta, S. K. Pati, and 



C. N. R. Rao, [Angew. Chem. Int. Edit., 49, 4059 (20To] 
i'^ G. Eda, H. Yamaguchi, D. Voiry, T. Fujita, M. Chen, and M. Chhowalla, [Nano Lett., 11, 5111 
(2011) . 



14 



20 



T. Korn, S. Heydrich, M. Hirnier, J. Schmutzler, and C. Schuller, [Appl. Phys. Lett., 99, 102109 



(2011) 



21 S. Lebegue and O. Eriksson, |Phys. Rev. B, 79, 115409 (2009) 

22 A. Splendiani, L. Sun, Y. Zhang, T. Li, J. Kim, C.-Y. Chim, G. Galli, and F. Wang, f^^ 
Lett., 10, 1271 (2010). 



23 J. K. Ellis, M. J. Lucero, and G. E. Scuseria, |Appl. Phys. Lett., 99, 261908 (20ll) 

2'* D. S. Inosov, V. B. Zabolotnyy, D. V. Evtushinsky, A. A. Kordyuk, B. Biichner, R. Follath, 



H. Berger, and S. V. Borisenko, [ New J. Phys., 10, 125027 (20 08)1 
25 T. Kidd, D. Klein, T. Rash, and L. Strauss, |Appl. Surf. Sci., 257, 3812 (2011 



26 A. Colev, C. Gherman, V. Mirovitskii, L. Kulyuk, and E. Fortin, [TLuniin., 129, 1945 (2009) 
2'^ B. Radisavljevic, A. Radenovic, J. Brivio, V. Giacometti, and A. Kis, Nat. Nanotechnol., 6, 



147 (2010) 



28 L. Liu, S. Kumar, Y. Ouyang, and J. Guo, [IEEE T. Electron Dev., 58, 3042 (20lT) . 
2^ Y. Yoon, K. Ganapathi, and S. Salahuddin, Nano Lett., 11, 3768 (2011). 



30 B. Radisavljevic, M. B. Whitwick, and A. Kis, |ACS Nano, 5, 9934 (20lT) . 

31 Z. H. Ni, T. Yu, Y. H. Lu, Y. Y. Wang, Y. P. Feng, and Z. X. Shen, |ACS Nano, 2, 2301 (2008") 



32 E. McCann, phys. Rev. B, 74, 161403 (2006) . 

33 E. V. Castro, K. S. Novoselov, S. V. Morozov, N. M. R. Peres, J. M. B. Lopes dos Santos, 



J. Nilsson, F. Guinea, A. K. Geim, and A. H. C. Neto, |Phys. Rev. Lett., 99, 216802 (2007) . 

34 T. Ohta, A. Bostwick, T. Seyller, K. Horn, and E. Rotenberg, [Science, 313, 951 (2006") , 

35 Y. Zhang, T.-T. Tang, C. Girit, Z. Hao, M. C. Martin, A. Zettl, M. F. Crommie, Y. R. Shen, 



and F. Wang, [Nature, 459, 820 (2009) . 
36 F. Xia, D. B. Farmer, Y.-m. Lin, and P. Avouris, [Nano Lett., 10, 715 (2010] , 
3^^ A. Ramasubramaniam, D. Naveh, and E. Towe, [Phys. Rev. B, 84, 205325 (20111 

38 E. Scalise, M. Houssa, G. Pourtois, V. Afanas'ev, and A. Stesmans, [Nano Res., 5, 43 (2012) 

39 M. Dave, R. Vaidya, S. G. Patel, and A. R. Jani, Buh. Mater. Sci., 27, 213 (2004). 

40 T. Boker, R. Severin, A. Miiller, C. Janowitz, R. Manzke, D. Vofi, P. Kriiger, A. Mazur, and 



J. Pohmann, [Phys. Rev. B, 64, 235305 (2001) 



41 K. F. Mak, C. Lee, J. Hone, J. Shan, and T. F. Heinz, [Phys. Rev. Lett., 105, 136805 (2010^ 

42 W. Schutte, J. D. Boer, and F. Jehinek, |J. Solid State Chem., 70, 207 (1987) 

43 P. E. Blochl, [Phys. Rev. B, 50, 17953 (1994^ . 



15 



^^ G. Kresse and D. Joubert, phys. Rev. B, 59, 1758 (1999) . 

^^ J. P. Perdew, K. Burke, and M. Ernzerhof, phys. Rev. Lett., 77, 3865 (1996) 

'^^ G. Kresse and J. Hafner, |Phys. Rev. B, 47, 558 (1993) 

47 



S. Grimme, [Journal of Computational Chemistry, 27, 1787 (20061 , ISSN 1096-987X. 
'^^ L. Wei, C. Jun-fang, H. Qinyu, and W. Teng, piiysica B, 405, 2498 (20T0] . 
"^9 W. Tang, E. Sanville, and G. Henkelman, | jrPhys.: Condens. Matter, 21, 084204 (2009) 
50 E. Sanville, S. D. Kenny, R. Smith, and G. Henkelman, |J. Comput. Chem., 28, 899 (2007) 



5^ G. Henkelman, A. Arnaldsson, and H. Jonsson, Comp. Mater. Sci., 36, 354 (2006) 

52 J. p. Perdew and A. Zunger, phys. Rev. B, 23, 5048 (198T 

53 J. Heyd, G. E. Scuseria, and M. Ernzerhof, |J. Chem. Phys., 118, 8207 (2003) 



54 J. Heyd, G. E. Scuseria, and M. Ernzerhof, |J. Chem. Phys., 124, 219906 (2006) . 

55 Y.-i. Matsushita, K. Nakamura, and A. Oshiyama, phys. Rev. B, 84, 075205 (201i; 

56 L. Hedin, phys. Rev., 139, A796 (1965" 



16 



